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Since 2006, EPOS hadronic interaction model is being used for very high energy cosmic ray analysis. Designed for 
minimum bias particle physics and used for having a precise description of SPS and RHIC heavy ion collisions, 
EPOS brought more detailed description of hadronic interactions in air shower development. Thanks to this 
model it was possible to understand why there were less muons in air shower simulations than observed in real 
data. With the start of the LHC era, a better description of hard processes and collective effects is needed to 
deeply understand the incoming data. We will describe the basic physics in EPOS and the new developments 
and constraints which are taken into account in EPOS 2. 



1. INTRODUCTION 

Air shower simulations are a very powerful tool to 
interpret ground based cosmic ray experiments. How- 
ever, most simulations are still based on hadronic in- 
teraction models more than 10 years old. Much has 
been learned since, in particular due to new data avail- 
able from the SPS and RHIC accelerators. 

In this paper, we discuss the new development of the 
EPOS model, the latter one being a hadronic interac- 
tion model, which does very well compared to RHIC 
data a, and also other particle physfcs experiments 
] (especially SPS experiments at CERN). Used for air 
shower analysis since 2006, the last version 1.99 re- 
leased in 2009 gives very interesting results in terms of 
mass composition . Due to the constraints of parti- 
cle physics, air shower simulations using EPOS present 
a larger number of muons at ground It allows 
for the first time to reproduce both the muon num- 
[ ber and the elongation rate using a reasonable aver- 
age mass for all energies between the knee and the 
Greisen-Zatsepin-Kuzmin (GZK) cut-off. 

On the other hand, the new measurements at LHC 
give us the opportunity to test the model more deeply. 
After a general introduction of the EPOS model, we 
will explain in this paper the new developments done 
in EPOS 2. 



2. EPOS MODEL 

One may consider the simple parton model to be the 
basis of high energy hadron-hadron interaction mod- 
els, which can be seen as an exchange of a "parton 
ladder" between the two hadrons. 

In EPOS, the term "parton ladder" is actually meant 
to contain two parts the hard one, as discussed 
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Figure 1: Elementary parton-parton scattering: the hard 
scattering in the middle is preceded by parton emissions 
attached to remnants. The remnants are an important 
source of particle production even at RHIC energies. 



above, and a soft one, which is a purely phenomeno- 
logical object, parameterized in Regge pole fashion. 

In additions to the parton ladder, there is another 
source of particle production: the two off-shell rem- 
nants, see fig. [H We showed in ref. [5:] that this 
"three object picture" can solve the "multi-strange 
baryon problem" of conventional high energy models, 
see ref. [1]. 

Hence EPOS is a consistent quantum mechanical 
multiple scattering approach based on partons and 
strings 4], where cross sections and the particle pro- 
duction are calculated consistently, taking into ac- 
count energy conservation in both cases (unlike other 
models where energy conservation is not considered 
for cross section calculations ,7]). Nuclear effects re- 
lated to Cronin transverse momentum broadening, 
parton saturation, and screening have been introduced 
into EPOS [8]. Furthermore, high density effects lead- 
ing to collective behavior in heavy ion collisions are 
also taken into account Q. 
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Thanks to a Monte Carlo, first the coUision config- 
uration is determined: i.e. the number of each type of 
Pomerons exchanged between the projectile and tar- 
get is fixed and the initial energy is shared between 
the Pomerons and the two remnants. Then particle 
production is accounted from two kinds of sources, 
remnant decay and cut Pomeron. A Pomeron may be 
regarded as a two-layer (soft) parton ladder attached 
to projectile and target remnants through its two legs. 
Each leg is a color singlet, of type qq , qqq or qqq from 
the sea, and then each cut Pomeron is regarded as two 
strings, Cf. fig. [2] a) and b). 
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Figure 2: a) Eacli cut Pomeron is regarded as two 
strings b). c) The most simple and frequent collision 
configuration has two remnants and only one cut 
Pomeron represented by two q — q strings, d) One of the 
q string ends can be replaced by a qq string end. e) With 
the same probability, one of the q string ends can be 
replaced by a qq string end. 

It is a natural idea to take quarks and anti-quarks 
from the sea as string ends for soft Pomeron in EPOS, 
because an arbitrary number of Pomerons may be in- 
volved. 

Thus, besides the three valence quarks, each rem- 
nant has additionally quarks and anti-quarks to com- 
pensate the flavors of the string ends, as shown in 
fig- m c). According to its number of quarks and 
anti-quarks, to the phase space, and to an excita- 
tion probability, a remnant decays into mesons and/or 
(anti)baryons Furthermore, this process leads to 
a baryon stopping phenomenon in which the baryon 
number can be transfered from the remnant to the 
string ends (for instance in [2] d) , depending on the 
process, the 3q -f 3q can be seen as 3 mesons or a 
baryon-antibaryon pair) . 

In case of meson projectile, this kind of diquark pair 
production at the string ends leads to an increase of 
the (anti)baryon production in the forward production 
in agreement with low energy pion-nucleus data IC [ 
as shown fig. [31 Comparing to QGSJEtII model |ll [ 
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Figure 3: Model comparison: longitudinal momentum 
distributions of pion carbon collisions at 100 GeV from 
EPOS with (full) or without (dashed-dotted) sting-end 
diquarks and QGSJETII (dashed) compared to data [lo| . 



which do not have diquark as string ends or using only 
qq as string end in EPOS, we can clearly see that this 
process is needed to reproduce experimental data. As 
a consequence it is part of the larger number of muons 
in air shower simulations with EPOS. 

Energy momentum sharing and remnant treatment 
are the key points of the model concerning air shower 
simulations because they directly influence the multi- 
plicity and the inelasticity of the model. 



3. NEW DEVELOPMENTS 

With the start of the LHC era, it is now possi- 
ble to develop some particular physics point of the 
EPOS model, which will be possible to test with newly 
available data. 



3.1. Hydrodynamical Evolution 

In EPOS 2, a new tool has been developed for treat 



12[ for details 



ing the hydrodynamical evolution (see 
and tests with AuAu data). 

As we saw, in case of elementary reactions like pro- 
ton proton scattering (at moderately relativistic ener- 
gies), hadron production is realized via string break- 
ing, such that string fragments are identified with 
hadrons. When it comes to heavy ion collisions or 
very high energy proton-proton scattering, the proce- 
dure has to be modified, since the density of strings 
will be so high that they cannot possibly decay inde- 
pendently. For technical reasons, we split each string 
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into a sequence of string segments, at a given proper- 
time To, corresponding to widtlis 6a and (5/3 in the 
string parameter space. One distinguishes between 
string segments in dense areas (more than some crit- 
ical density po of segments per unit volume), from 
those in low density areas. The high density areas are 
referred to as core, the low density areas as corona 
. String segments with large transverse momentum 
(close to a kink) are excluded from the core. Based 
on the four-momenta of infinitesimal string segments, 
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one computes the energy-momentum tensor and con- 
served currents. The corresponding energy density 
e(ro,x) and the flow velocity v{to,x) serve as initial 
conditions for the subsequent hydrodynamic evolu- 
tions, which is characterized by : 

• consideration of the possibility to have a (mod- 
erate) initial collective transverse flow; 

• event- by- event procedure, taking into the ac- 
count the highly irregular space structure of sin- 
gle events, being experimentally visible via so- 
called ridge structures in two-particle correla- 
tions; 

• use of an efficient code for solving the hydrody- 
namic equations in 3-1-1 dimensions, including 
the conservation of baryon number, strangeness, 
and electric charge; 

• employment of a realistic equation-of-state, 
compatible with lattice gauge results - with a 
cross-over transition from the hadronic to the 
plasma phase; 

• use of a complete hadron resonance table, mak- 
ing our calculations compatible with the results 
from statistical models; 



orders too. Since this object is very difficult to cal- 
culate explicitly at all orders including proper energy 
conservation, we use a simple parameterization hav- 
ing the same form as the usual Pomeron (see 8]) and 
fixing the parameters using the proton-proton single 
difFractive cross section and the energy spectrum of 
leading protons. 

3.3. Structure Function 

In order to test hard scattering at LHC (or Teva- 
tron) using EPOS, it is important to be able to repro- 
duce the structure function F2 as measured thanks to 
the HERA electron-proton collider. 
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Figure 4: "Half-" Pomeron corresponding to the parton 
distribution function tested via the structure function F2. 

In EPOS, F2 corresponds to a "half" hard or semi- 
hard Pomeron as illustrated in fig.lH The cross section 
of a quark with a given as a function of the mo- 
mentum fraction x can be calculated using the convo- 
lution of a soft pre-evolution with the DGLAP equa- 
tions for the perturbative development of the parton 
ladder Taking into account the effective correc- 
tions due to higher order diagrams in the connection 
between the Pomeron and the nucleon (and fixed with 
proton- proton cross section), we obtain the structure 
function as shown in fig. [Sj 



• hadronic cascade procedure after hadronization 
from the thermal system at an early stage. 

3.2. Diffraction 

In order to produce hard diffractive events as mea- 
sured at Tevatron and most likely at LHC in the near 
future, it was necessary to improve the way of produc- 
ing difFractive events in EPOS 2. To get a consistent 
description of low mass and high mass diffraction, an 
effective diffractive Pomeron is used in addition to the 
soft and semi-hard Pomeron. This effective object rep- 
resents in reality all type of diagrams (soft and hard 
component) which are not connected to at least one 
of the remnant, leaving the latter intact. The triple 
Pomeron would be part of this object, but all higher 



3.4. Baryon production 

In EPOS 2, not only the soft Pomcrons can have a di- 
quark as string-end, but the semi-hard Pomerons are 
now treated the same way. It will increase the produc- 
tion of (anti-)baryons in the forward region compared 
to the previous version of EPOS. This may be seen in 
the number of muons produced in air shower simula- 
tions. 

Furthermore, to take into account the property of 
different jet types, the string tension used for the 
string fragmentation now depends on the initial par- 
tons. Jets coming from quarks or from gluons will not 
produce the same ratio of (anti-)baryons or strange 
particles over pious, and it will change the energy evo- 
lution of such ratios. 
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the better calculation of the structure function and 
the more accurate treatment of the hydrodynamical 
phase, provide much better results in comparison with 
LHC data than with EPOS 1.99. We show the pseudo- 
rapidity distribution at 7 TeV in fig. [6] as an example. 
EPOS 2 results are compared to epos 1.99 and other 
models used for air shower simulations QGSJEtII and 
Sibyll 2.1 [H and ALICE data [l^. 

More detailed results can be found in [l^, [l^ , es- 
pecially concerning the effect of the hydrodynamical 
evolution on correlations between secondary particles 
(Bose-Einstein correlation, ridge). 

For the moment EPOS 2 can only be used for mini- 
mum bias physics, but in the future it will be possible 
to select special class of hard events to study specific 
channels and underlying events (20l |. 
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Figure 5: The structure function F2 for different values 
of together with experimental data from HI [ij 
ZEUS flil and NMC 
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Figure 6: Pseudorapidity distribution in pp scattering at 
7 Te'V (INEL>0 trigger), compared to ALICE data 
(points) for different models : results from EPOS 2 (full) 
together with results from EPOS 1.99 (dashed), QGSJEtII 
(dotted) and Sibyll 2.1 (dashed-dotted). 
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